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Abstract 

The light-cone wave function of the nucleon is calculated in the limit N c — > oo in 
the quark-soliton model inspired by the theory of the instanton vacuum of QCD. The 
technique of the finite time evolution operator is used in order to derive expressions 
for all components of the Fock vector describing the nucleon in the infinite momentum 
frame. It is shown that nucleon wave function for large N c can be expressed in terms of the 
wave function of the discrete level in the self-consistent meson field and light cone wave 
functions of 1,2, etc mesons. The 3-quark components of the nucleon and A-resonance 
are estimated. Wave function of the nucleon appears to be positive in the whole region of 
x and it has rather small asymmetry. It differs strongly both from Chernyak-Zhitnitsky 
wave function and the asymptotic one. Large momentum transfer asymptotic of the 
electromagnetic and axial form factors is discussed. 



1 Introduction 

Light cone wave functions of hadrons were introduced in hadron physics many years ago 
[H [21 EE El 13 They contain virtually the full information which is necessary to describe 
hadron properties at high energies. The use of hadron wave functions in the context in 
QCD relies on the concept of factorization. Processes with hadrons are divided into 2 
parts: i) the hard pro cess- dependent parts are calculated according to perturbative QCD 
and ii) the soft process-independent part is usually encoded in soft functions, parton 
distributions, fragmentation functions, etc. Usually, a formal definition of the soft part 
can be formulated in terms of definite quark and gluon hadron matrix elements. Their 
logarithmic scale dependence is well understood in terms of corresponding evolution equa- 
tions. 



In principle, any soft part can be expressed in terms of the hadron wave function. For 
example, ordinary parton distributions is a sum of wave functions squared corresponding 
to a different number of quarks or gluons integrated over all momenta of partons except 
one. The analogous definition in terms of wave functions can be done for the fragmentation 
functions as well. 

Another type of hadronic matrix elements is involved in the description of the elastic 
form factors or transition form factors. There one deals with the matrix elements of quark 
currents between hadron states which have very different initial and final momenta. The 
large difference of momenta in exclusive reactions effectively separates out the component 
of wave function with the minimal number of constituents, so gives an access to the 
simplest structures in hadron jTJ 13 E3 HI E] • 

An intermediate situation between these two types of reactions occurs in processes like 
deeply virtual Compton scattering and hard meson production. In this case, the hadron 
operators involved are bilocal and their matrix elements are off-forward. The soft part of 
such processes is described by so-called generalized distributions. Recently it was shown 
that generalized distributions can also be easily presented in terms of light cone wave 
functions 0. 

Unfortunately up to now light-cone wave functions of baryons (as well as mesons) 
in the low normalization point cannot be determined from the first principles of QCD. 
Perturbative theory is able to predict only so-called asymptotic wave functions which are 
normalized at arbitrary high normalization point. 

Up to the present moment the light-cone nucleon wave functions were calculated only 
in the framework of QCD sum rules [7j. Of course, this can be considered only as a crude 
estimate. For this reason self-consistent models of the nucleon (more or less motivated 
by QCD) become highly desirable. Unfortunately, to the best of our knowledge, no 
calculation of the nucleon wave function was made in any self-consistent relativistic field- 
theoretical model. 

In this paper we attempt to calculate light-cone wave functions at a low normalization 
point in the limit of large number of colours, N c — > oo. Even though in reality N c = 3, 
the academic limit of large N c is known to be a useful guideline. It is a general QCD 
theorem that at large N c the nucleon is heavy and can be viewed as a classical soliton jHj. 

An example of the dynamical realization of this idea is given by the Skyrme model jHj- 
However, the Skyrme model is based on an unrealistic effective chiral Lagrangian. A far 
more realistic effective chiral Lagrangian is given by the functional integral over quarks 
in the background pion field [TU] : 

exp (iS cS [n(x)}) = J DtfjD^exp (i J d 4 xtp(i$ - MW 5 )^ , 

U = exp (m a (x)T a ) , f/ 75 = exp (i^ a (x)r a -f 5 ) = ~^-U + ^—^U ] . (1) 

Here ip is the quark field, M is the effective quark mass which is due to the spontaneous 
breakdown of chiral symmetry (generally speaking, it is momentum-dependent) and U 
is the SU(2) chiral pion field. The NJL-type Lagrangian Eq. (P) has been derived from 
the instanton model of the QCD vacuum ^T] FEZ] , which provides a natural mechanism 
of chiral symmetry breaking and enables one to express the dynamical mass M and the 
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ultraviolet cutoff A intrinsic to eq. (JTJ) through the Aqcd parameter. Proposed on the 
basis of Lagrangian of eq. (JTJ) the chiral quark-soliton model |T2J ECU El UH] describes 
properties of baryons far better than the Skyrme model. For the recent status of the 
chiral quark-soliton model see reviews [HUE!- 

Recently in the framework of quark-soliton the parton distributions [T%1 ITU] and the 
generalized distributions [201 have been computed at a low normalization point. In this 
paper we complete the program of the investigation of high-energy properties of baryons 
in the quark-soliton model and present the framework for calculation of the light-cone 
wave of the nucleon. 

We use the evolution operator technique in order to present wave functions of the 
nucleon state moving with the speed V. In the limit V — > 1 the corresponding wave 
function tends to the light-cone wave function in question. This method is convenient as 
the evolution operator can immediately be expressed as a functional integral with definite 
boundary conditions. 

At large N c the functional integral representing the evolution operator can be evaluated 
using the saddle-point method. The saddle-point of the effective Lagrangian in eq. ((H) 
in the sector with unity baryon charge corresponds to the nucleon-soliton The 
stationary solutions of saddle-point equations correspond to the nucleon in the rest frame. 
Owing to relativistic invariance of the equations of motion there is also infinite number of 
other solutions which describe the moving solitons. We need the solution which describes 
the nucleon in the infinite momentum frame, and we have to extract the corresponding 
quark- ant iquark wave function. In the N c — > oo limit it can be viewed as a product of 
the quark states in the time-dependent self-consistent pion field. 

We show that the nucleon wave function at large N c is a product of N c one-quark wave 
functions of the valence quarks and the coherent exponential of quark- ant iquark pairs 
corresponding to the sea quarks. As it should be in this limit, the nucleon wave function 
is completely factorized in colours. One-quark wave functions receive contribution both 
from the wave function of the discrete level in the mean pion field and the sea quarks. 
The quark- ant iquark pair wave function can be expressed in terms of so-called Feynman 
Green function at finite time. 

This structure of the light-cone nucleon wave function is rather general. Indeed, as was 
already said, nucleon is a soliton of an effective meson Lagrangian at large N c . Hence its 
quark wave function is a product of states in the external field of all mesons. The model 
with Lagrangian of eq. (JTJ) (which is, in turn, based on the theory of the instanton vacuum) 
is specific in two respects. First, as the size of the nucleon (~ 1/M) is parametrically 
large as compared to the ordinary hadron scale (~ 1/p, p being the size of the instanton), 
only the lightest degrees of freedom (i.e. pions and constituent quarks) are important. 
Second, the number of gluons is suppressed in the instanton vacuum by the parameter 
(Mp) 2 <C 1 |TSJ E| ; so gluons in this model do not participate in the formation of the 
nucleon wave function. 

All what one needs to know in order to calculate the nucleon wave function in the 
model is the wave function of the discrete level (the solution of the Dirac equation in the 
external field) and the quark-antiquark pair wave function. Expanding the latter wave 
function in 7r-meson field, it is possible to present the wave function as convolution of the 
pion mean field with the quark-antiquark light-cone wave function of 1,2,. . . pions. This 
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is to be expected, as the pions are the only agents inducing the interaction in the model. 

In fact, light-cone wave functions of one and two pions were already considered in the 
instanton vacuum (see Refs. [2Sj and [21], correspondingly) and appeared to describe the 
data correctly. We discuss the connection with meson wave functions in the Section 3. 

We formulate the scheme for calculation of the nucleon-soliton light-cone wave function 
in Section 2 and in Section 3 we consider the wave functions in the infinite momentum 
frame. The three-quark components of the nucleon wave functions (so-called distribution 
amplitude) is calculated in Section 4. They appear to be almost symmetric (antisymmetric 
part of the nucleon wave function is numerically small) smooth functions, which are far 
from both the asymptotical wave function and the wave function of Chernyak-Zhitnicky 



We also discuss physical observables (asymptotics of electromagnetic and axial form 
factors) and evolution of the model wave function to the high Q 2 . Asymptotics of the 
form factors cannot be calculated explicitly as corresponding integrals depend strongly 
on the region of small x where the model is not valid. However we are able to calculate 
ratios of the form factors which appear to be rather close to the experimental data. 

2 Soliton wave function in field theory 
at large N c 

In principle, calculation of the wave function of a given state in terms of quarks and 
antiquarks is straightforward in the quantum field theory. However, usually this task is 
too complicated. In fact, we know only one model field theory, namely, the Schwinger 
model where the program of calculating of all wave functions was completed [21]. The 
most direct way to obtain wave functions of any state is to calculate the evolution operator 
S(T) for the given theory and present it as a sum: 



Here \n) is vector of certain state — the eigenfunction of Hamiltonian H . As it is well- 
known, the evolution operator of eq. (J2J) can be expressed as a functional integral at finite 
time with definite boundary conditions, namely 1 

1 We immediately write down the functional integral in the model of eq. Q but, in fact, we have to start 
from the full QCD. The corresponding functional integral is the integral over both the quark and gluon 
field. In the model of the instanton vacuum this integral is saturated by instanton and antiinstanton 
gluon configuration. One can then integrate out gluons in this approximation. As it was shown in 
Refs. ^3 El j one obtains the low-energy effective Lagrangian eq. JQ). As the number of gluons in the 
nucleon in the instanton model is parametrically suppressed, we should not put any boundary conditions 
on gluon field. For this reason the derivation of the evolution operator in the low-energy limit repeats 
literally the derivation of the low-energy effective lagrangian that was done in these papers. If one is 
still interested in gluon components of the wave function, they should be traced directly from the general 
expression for the full QCD evolution operator to which we have to apply the same approximations which 
lead to the effective Lagrangian of eq. This is straightforward in the instanton vacuum model. Gluon 
components are of the order of {M p) 2 <C 1 where p is the instanton size. The same parameter was used 
in the derivation of effective lagrangian eq. Q 



type [29 E2| • 




(2) 
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S[T] = / 7^(x) / Di,{x) / Dtt(x) exp i dt £ e// , (3) 

where £ e // is the effective Lagrangian of eq. (0). The operators a ± (p), b ± (p) are annihilation- 
creation operators for quarks-antiquarks which are defined through the expansion of the 
field ip(x) into positive- and negative-frequency parts: 

M*) = /t^Ei/^ U x \p)u (x \p)e^ + b^ + {p)v^\p)e-^) (4) 
J {in) 6 x V Po 

(here A is quark polarization; u^ x \p), v^ x '{p) are free plane wave spinors). 

We do not impose any boundary conditions on the 7r-meson field. This field appears in 
the derivation of the effective Lagrangian as a result of bosonization JU|, and it should not 
be considered as an elementary one. In fact, it is impossible to consider the light-cone wave 
function in terms of both quarks and 7r-mesons — it would be a kind of double-counting. 

In the representation of eq. (jHJ the quark-antiquark creation operators a + ,b + anti- 
commute with annihilation operators a, b, as they are related to different times and 
T. One has to calculate the functional integral of eq. (J3J) at the finite time T as a func- 
tional of operators a ± , b ± . Expanding in the exponentials exp(—iE n T) and factorizing 
creation operators from annihilation ones, we can find eigenfunctions of all states in terms 
of quarks-antiquarks (see, e.g. [2T]). 

In the large iV c -limit the functional integral over the 7r-meson field should be calculated 
in the saddle-point approximation. This procedure can be formulated as follows. 

Let us first integrate over the quark field ip(x). We divide fermion field into two parts, 
ip(x) and xi x )'- 

if>(x)=i}>(x)+x(x), ${x) = $(x) +x(x). (5) 

Let us require that ip(x) obeys the Dirac equation in the external 7r-meson field with 
corresponding boundary conditions: 

(id- Me™ {xh5 )i> = 0, 
& + \2,t = 0) = A(x) = J ^e^^^a x {p)u\p), 

&-\2,t = T) = B+{x) = J ^e-^y^E^(p)^(p), (6) 

The field x( x ) obeys zero boundary conditions. 

A solution of the Dirac eq. © can be expressed in terms of the finite time Feynman 
Green function G^ T \x, t\y, 0) 

d 3 p 



${x,t) =Jd 3 yJ ^y/§E [G {T \x,t\y,0)u x (p)a x (p)e^y + 

+G^ T \x,t\y,T)v\p)bi(p)e-^y] , (7) 
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with an analogous expression for ip(x, t). 

The finite time Green function G^ T \x, t\x', t') is the solution of the Dirac equation 
in the external field which has vanishing negative-frequency part at t — and vanishing 
positive-frequency part at t = T. At T — > oo finite time Green function reduces to the 
usual Feynman Green function. Feynman propagators for free particles appear to obey 
automatically the required boundary conditions at any T. Also it is easy to verify that the 
finite time Green function in the external field can be constructed as a sum of diagrams 
with free Feynman propagator interacting with external field. The only difference with 
the usual Feynman Green function is that all integrals in intermediate points should be 
carried only in the finite time interval < t < T (for Feynman Green function one has to 
integrate over the whole interval — oo < t < oo). 

Let us substitute the expansion of eq. J5J) into the functional integral of eq. (jHJ). One 
can see that variables if) and x( x ) are completely separated from each other. As a result 
we arrive at the following expression for the evolution operator (cf [2T]): 

5[r,a ± ,6^] = jD7i(x)Det {T) [7c]S ext {7c,T] (8) 

(we used here eq. (J7J) for if)). Here S ex t[n, T] is the evolution operator in the given external 
field: 

S ext [7i,T] = expyd 3 xd 3 y[B(x) 1 G m (xO,yT)B + (y) + 

+B{x) 1 °G iT \xe, yO)A{y) + A + (x)^ G {T \xT , yO)A{y) + 

+A + (x) 1 °G^ T \xT - e, yT)B + {y)] } , (9) 

where e —>■ +0); the quantity Det^[7r] is the finite-time determinant in the external field. 
It is the Gaussian functional integral over the fermion field x( x ) with zero boundary 
conditions. For this reason it does not depend on the operators a , b ± being only the 
functional of pion field tt(x). It can also be presented in terms of the finite-time Green 
function: 

rM pT r 

dM dt d 3 xTr [G {T) (xt, xt)iU 15 (x, t) y 



Det^ 



7rj = exp 



o Jo 



(10) 



Thus the evolution operator in the external field is the coherent exponential of the 
creation-annihilation operators. The nucleon is the lowest possible state in the sector 
with baryon charge B — 1. One can obtain its wave function by applying the evolution 
operator to any colourless state of N c quarks and taking the T — > oo limit. For example, 
the nucleon wave function <3>tv can be obtained from the state of free quarks: 



N c 

c N e- iM » T <S> N (a + ,b + )=\im 5(T) J] <(Pi) l«o> 



i=l 

N c 



~ jD^et^}\{GT\pi^h,T)at{k)eM^{p) G (T) (T,p, T,p') b + (p')] (11) 

< 

(we write the formula in a bit symbolical form, at the moment it is sufficient for our 
purposes). Here M.n is the nucleon mass and cn is the overlap between the initial 
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state of N c quarks and nucleon wave function. The functional integral of eq. (|11|) should 
be calculated in the saddle-point approximation. One has to find the pion field which 
extremizes the integrand. At large T important factors are 2 

Det (T) [vr] ~ exp(-iE field [7r]T), G? \ Pi , 0, k h T) ~ exp(-iE level [ir}T) , (12) 

where Ef^i^} is the energy of the Dirac continuum with the given pion field (proportional 
to N c ) and Ei eve i [ir] is the energy of the (possible) discrete level for the quark in this field. 
In order to find the saddle-point one has to minimize the sum: 

S[tt] = E field [7c] + N c E level [n] (13) 

in the presence of the pion field. It is exactly the condition which was used in constructing 
the nucleon-soliton in Refs ^21 EH E31 E3- Both contributions to the total energy are of 
the order of N c . The first is the full effective chiral Lagrangian (ECL) in the low-energy 
effective theory of eq. ([TJ calculated for a given pion field n. 

The minimum of eq. (|13|) is achieved at some stationary pion field and corresponds to 
the nucleon at rest. The value of the energy in the minimum is the nucleon mass: 

M N = 7nin£[n]~0(N B ). (14) 

This minimum was found in Ref. ^3 EI| and it corresponds to the hedgehog symmetry 
of the pion field: 

n a (x) = n a P(r), n = -, (15) 

r 

where the profile function P(r) is to be calculated numerically. 

We can obtain the wave function of nucleon in the leading order of N c , if we substitute 
the saddle-point field tt(x) of eq. (|15|) into eq. ([lip. In the higher orders, one has to express 
the general pion field n(x,t) = 7t + ■ n ( i uant anc i then perform the Gaussian integration in 
TT guant which appears to be ir guant ~ 0(l/\JW c ). Thus it is possible to formulate the 
systematic perturbation theory in 1/N C . However, in this paper, we restrict ourselves to 
the leading order. 

Let us stress that the procedure of the calculation of the nucleon wave function which 
we have formulated, is rather general. Indeed, it is a general QCD theorem that at large 
N c , the nucleon is the soliton of some effective meson Lagrangian. Thus its wave function 
can always be presented in the form of eq. (jlljl where Green functions should be found 
in the self-consistent field of all mesons entering this effective Lagrangian. Of course, 
the exact low-energy meson Lagrangian is unknown. In the present work, we use the 
instanton vacuum model in order to fix this low-energy Lagrangian. 

Let us also rewrite nucleon wave function in a different form. As explained above, 
nucleon can be described as N c valence quarks + Dirac continuum in the self-consistent 
external field. It is clear from eq. (fTTj) that wave function of the Dirac continuum (i.e. 

2 Operator exponential in eq. 1)11(1 does not contribute to the saddle-point eq. 1)13(1 as the Green function 
C7^ T '(T,p, T,p') does not contain exponential with the phase proportional to the time T. Let us note 
also that minimizing of eq. 1(1311 gives the saddle-point for pion field only at times which are far from the 
end-points of the interval (0,T). In fact this is enough for the calculation of wave function. 
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the state with all states with negative-energies occupied) is the coherent exponential of 
the quark- ant iquark pairs: 



| ft) = exp 



= exp 



]T J d 3 xd 3 y A + {x) 7°G (T) (xT - e, yT) B + (y) 
d 3 pid 3 p2 



.colour 



|fio> 



E 

.colour 



(2tt)* 



(16) 



where |Oo) is the vacuum of quarks-antiquarks. The function Al ' A2 (pi,p2) can be called 
the wave function of the quark-antiquark pair. We see that it is expressed in terms of the 
finite time Green function at equal times. 

The nucleon itself is a state with one more level occupied by valence quarks, namely, 
a discrete level (with positive energy) which appears in the self-consistent external field. 
We can obtain the nucleon wave function by applying the operator which fills this discrete 
level to eq. (jIBj): 

®n= II [ d 3 x^ + (x)f lev (x)\n), (17) 



colour 



where f\ ev (x) is the wave function of the discrete level (solution of Dirac equation in the 
external ir- meson field). In order to express $jv in terms of quark-antiquark we have 
to expand ^-operators according to eq. (J3J and commute them with the exponential of 
eq. ()16jl. As a result we get the following expression for the nucleon wave function (A is 
the normalizing constant): 



N 



x exp 



E 

.colour 



=a n 

colour 

d 3 pid 3 p 2 
(2vr) 3 



d 3 p 



a$(p)x 



(2tt) 3 

aUPi) e Al ' A2 (Pi,P 2 ) b+M) 



Here F x (p) is the one-quark wave function. It is a sum of two contributions: 



F\p) 



u*\p) fUP) (2?r) 3 5 (3) (p-f)- 



d 3 p' M 
(2tt) 3 V 

-e^'(pj) v* x 'a/) M-f) 



(19) 



The first contribution in eq. (|T9|) is that of valence quarks, while the second term can be 
called the contribution of the sea quarks to one-quark wave function. 

Let us point out the fact that the nucleon wave function is completely factorized in 
colours. In fact, it is a general theorem in the strict N c —>■ oo limit. 



3 Wave function of the nucleon in the IMF 

It is well-known that the wave function in the rest frame has not too much physical sense 
jS]. The physical meaning can be ascribed only to the wave function of the fastly moving 
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nucleoli (nucleoli in the infinite momentum frame). Contrary to the wave function in the 
rest frame it can be accessed by measurements. 

The common approach to the hadron wave functions is the formalism based on the 
light-cone quantization (see, e.g., review j22] )• This formalism has many obvious advan- 
tages as compared to the approach based on Schrodinger wave functions. From the other 
side, the Schrodinger equal-time wave function can be defined in any frame, not necessar- 
ily on the light-cone. Also one can use well-trodden path which starts from the calculation 
of the evolution operator in the functional integral technique. This is why we use in this 
paper the last method. Light-cone wave function of the nucleon is, by definition, its wave 
function in the infinite momentum frame • 

The stationary saddle-point n a (x) corresponds to the nucleon-soliton at rest. However, 
as the effective chiral Lagrangian is relativistically invariant, we are guaranteed that 
there are infinitely many solutions of saddle-point equations of motion which describe 
the nucleon moving in some direction with a speed V. The corresponding pion field is 
time-dependent and can be obtained from the stationary field by Lorentz transformation 

^" = *(^f)- < 20 > 

In order to find all states of the Dirac continuum in the moving nucleon it is sufficient to 
solve the Dirac equation in the field of eq. (J2Uj) . In particular, the wave function of the 
valence level can also be obtained as Lorentz-transformation: 

<M*,f) = S[V)f lev (x) ^==f) exp (-i;^=2f j , (21) 

where e is the energy of the discrete level. Here S[V] is a matrix which transforms Lorentz 
indices 

S [V] = exp(ia 03 u), = -[7^,7^], tanh(^) = V. (22) 

The wave function of the discrete level in stationary field 7t(x) of eq. ()15|) is a mixture 
of two wave functions with orbital angular momentum L = 0, h(r), and L = 1, j(r) 



f*f(* - J_ ( ~ h (r)e af + ij(r) r -a^f . 



Here a is a spinor index and / is flavour index. The spherically symmetric functions f(r) 
and g(r) were found by numerical integration of the Dirac equation. 

Let us begin with the calculation of the contribution of the discrete level in the one- 
quark wave function (first term in eq. p9jl ) in the infinite momentum frame (V — > 1). 
We proceed in eq. (|21j) to momentum space and obtain from eq. (|19jk 

fleUk)= [ d 3 x e-^f«l(x). (24) 
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Let us divide the quark momentum into the longitudinal and transverse parts with respect 
to the total nucleon momentum P N = M. N V/ yl — V 2 : 

p = (zP N ,p±) . 

Also in the IMF it is convenient to use quark- ant iquark operators normalized in a different 
way than the usual ones: 



{a + (z 1 ,p 1± ),a(z 2 ,p 2± )} = 8{zi- z 2 ) (2tt) 2 5 (2) (p 1± - p 2 ±) 



(25) 



(We shall call a wave function the coefficient in front of quark operators normalized this 
way). We have from eq. 



Mm M 



Pn V z 



u(p)S[V]fUP)\ 



P3=zMn~£- 



(26) 



We make use of the relation 7 S[V"] = S 1 [V] , ~f and apply Lorentz transformation to the 
free spinor u{p): 



u {x \p)S- l [V]=u A (p), p 3 



P3 + Vu p 



Up + VP3 



(27) 



here u p = y/p 2 + M 2 . In the limit V — > 1, the momentum p 3 is large, p% ~ lu p ~ 
2zP^/ ' M.n- Wave functions of free fast-moving quarks become eigenfunctions of polar- 



ization operator 7°7 3 : 



—u^(l + tV), «W(p) = P x l^v^{l + 7 V), (2f 
m M N V m MN 



fzW(p) = P. 

V m MN 

where the quark- ant iquark spinors with spin up and down are: 



T 



1 



V2 







f° \ 















-1 


T 1 





'^"v/2 


1 


-1 







1 





\o j 




1 1 J 




I o y 




I 1 / 



(29) 



In fact, the matrix |(1 + 7°7 3 ) is a projector on the state with definite chirality. 

Finally we get for the contribution of valence quarks to the one-quark wave function 



FV(z,p x )=fUp±,*M N -e) 



(30) 



where ff ev {p) is Fourier transform of the level wave function. Using expression of eq. (|23j) 
we obtain: 



p2 



(W(p) -ph(p))<T3 + j(p)p±a± 



| P3=zM N -£ 



(31) 



where 



Hp) 



drr h(r)\ — 



2 sin pr 



it r 



j(p) = 
10 



drr j(r)\j — 



2 sin pr — pr cos pr 



pr^ 



We see that, as it should be, the dependence on the nucleon momentum P is cancelled 
in the final expression for the level contribution of eq. (|31|) to the nucleon wave function. 
Also it is clear that this contribution remains stable in the large-iV c limit in the main 
region of the quark momenta: Zi ~ 0(1/ N c ), p± ~ 0(1) (nucleon mass is proportional 
to the number of colours iV c , factor \/M.n is needed for the correct normalization of the 
one-quark wave function). 

Let us proceed now with the sea quark contribution to the wave functions. According to 
eq. (|19|) it can be expressed in terms of the quark-antiquark pair wave function A ' A (p, jf) 
which, in turn, is calculated through the finite-time Green function in the external pion 
field. 

The equal-time Green function is a sum of diagrams in the pion field. These diagrams 
are similar to the ordinary Feynman graphs except for the fact that the integration in the 
intermediate times in the graph is going only in the interval (0, T) (usually, one integrates 
over the whole interval (— oo, oo)). It is easy to show that these diagrams obey both the 
Dirac equation and necessary boundary conditions at t — and t = T. 

The diagrams represent the Green function as an expansion in the powers of the pion 
field. One can rearrange this series in such a way that the expansion is carried out in 
the powers of (U 15 — 1) instead. It can be shown that the new expansion is in increasing 
powers of the gradients of the pion field. 

Let us restrict ourselves by the first non-trivial term of this expansion. We call this 
approximation for the Green function an interpolation formula |T3*1 ITl] . It becomes exact 
in three limiting cases of the pion field: i) pion field is small ir(x) <C 1, ii) pion field 
is slowly varying function of coordinates with typical momenta k^ <C M, iii) pion field 
momenta are large k n 3> M. As a result, interpolation formula usually works rather well 
in the problems related to nucleon-soliton: its accuracy is typically around 10%. 3 

We introduce the Feynman Green function in the mixed p, t representation: 

rf^A ^ P 7°sign(t) - p7 / 

G(p,t) = — , u p = ^Jp 2 + M 2 . (32) 

The first-order correction to the pair wave function in the pion field in the IMF is equal 
to 

©Ai,A 2 (Pi,P2) = J d 3 k^^5 ^ -j=t== - ( Pl ) z - (p 2 ))j 5 (k ± -px± -p 2 ±) x 

x u^( Pl h g + ( Pl ) M < k \ g-(- P2 )v x *( P2 ), (33) 

signs ± label the sign of t in the free Green functions. Let us denote by z\ and z 2 the 
fractions of the nucleon momenta carried by the quark and antiquark correspondingly: 
(Vi,2)z = -21,2-P/v- The energy denominator in eq. (J3l!j) is small only if 1 > z\ 2 > 0. 



3 We checked the accuracy of this formula, for example, in the calculation of the nucleon structure 
function |18j confronting it to the exact calculation accounting for all solutions of the Dirac continuum. 
This problem is, of course, rather close to the calculation of wave function. The interpolation formula 
was accurate enough. 
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Expanding it in the powers of nucleon momentum and summing over spinor indices with 
help of eq. ((23), we obtain: 



©A^AaOi, z 2 ,p±i,p±a) = \M{z\ + z 2 )(cr 3 )x 1 \ 2 + {z 2 Px - ^iP 2 L )( ( ?" L )aiA 2 ) x 

M N M[Tt{pi +P2)]p lz =zxM N , P 2z=z 2 M N 



(34) 



(*i + z 2 ) Zl z 2 M 2 N + (M 2 + p\ 1 )z 2 + (M 2 + p\ 2 ) Zl 

where a, are Pauli matrices. 4 As it should be, the nucleon momentum P^ is cancelled 
out in this expression. We see again that pair wave function has a finite limit at large N c 
— in the same sense as the one-quark wave function. 

Collecting the terms in the perturbative expansion which correspond to the expansion 
of U 15 (x) in powers of n we can promote eq. (jSlj) to the interpolation formula. It has the 
following form: 

Q ( , MM N 

Q\ u x 2 {z 1 ,z 2 ,p± 1 ,p± 2 



{z 1 +z 2 )z 1 z 2 M 2 N + (M 2 +p 2 ±1 )z 2 + (M 2 +p 2 ±2 )z 1 
x {{Mfa + z 2 )(a 3 ) Xl \ 2 + (z 2 p±i - z l p ±2 )(a ± ) XlX ^j fi(pi +p 2 )- 
-i (M(z 2 - Zl )6 Xl \ 2 +i£ap(z 1 p 2 L a -z 2 p± a )(o-^) Xl \ 2 ) S(pi+p 2 )]} • (35) 



Here S(/c) is the Fourier transform of the scalar component of U l5 (x) and H(k) is its axial 
component. In the soliton mean field of eq. (|THj) they are equal to 

t(k) = J a?xe-^ (cos P(r) - 1) , fl(k) = J d 3 xe^ {nr) sin P{r), (36) 

The expression for the pair wave function of eq. has a clear physical meaning. The 
factor in front of the parenthesis in this expression is, in fact, a light-cone wave function 
of the pion in the leading order in iV c calculated in Ref. (2H]- 5 According to eq. (|55j) 
one has to convolute this wave function with a mean pion field in the nucleon in order to 
obtain quark- ant iquark pair wave function. This factorization is a natural consequence 
of the large N c limit and strictly corresponds to the soliton nature of the nucleon in this 

4 We normalize eq. I|34(l in such a way that the pair wave function is equal to 

dz 1 d 2 p 1± f dz 2 d 2 p 2 ± +(A ), , , m,+(A 2 )/ \ 



(2tt) 2 J (2tt) 2 
where quark operators are normalized to 

{a+^(z 1 , Pl± )a^(z 2 ,p 2± )} - S XlX2 5( Zl - z 2 )(2^) 2 ( 5« (p 1± - p 2± ) 

{b + ^\zi,pi±)b^(Z2,P2x)} = S XlX2 S( Zl - Z 2 )(2tt) 2 S^( Pi1 _ -P2A.) 



5 Pion light-cone wave function of Ref. |23|is obtained by substitution Z\ + z 2 = 1 (fractions should be 
measured relative to the momentum of the pion) and integrating over transverse momenta Pi±,p 2 ±, with 
the condition Pi±+p 2 j_ = 0. It corresponds to one definite (75) component of the full wave function. Also, 
in Ref. we took into account the dependence of the dynamical quark mass on the quark virtuality 
which is neglected here. 
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limit. Moreover, considering next terms in the expansion of the Green function one can 
prove this factorization also in all orders. For example, the second term is the product of 
the two-pion light-cone wave functions (calculated in our model in Ref. j21]) and the pion 
mean field squared. At last, according to eq. ()19j). the sea contribution to the one-quark 
wave function is a convolution of pair wave function with the antiquark contribution to 
the light-cone wave function of the discrete level. This formula is also guaranteed by the 
large N c limit. 

Projecting the conjugated wave function of the level onto antiquarks we obtain anal- 
ogously to eq. (j3Tj) 



F/* val (P±, z) = ~ VAW A) (1 + 7°7 3 )/iev(-p±, ~zM N - e) 



p 2 



(Psj(p) - ph(p)((Ti){ + j(p) (Px(v3){ + ip y 8\)] • (37) 



In order to obtain the sea contribution to the 1-quark wave function we have to convolute 
the latter with the pair wave function of eq. (j3*3j) : 



Again this contribution does not depend on N c in the same manner as the level wave 
function of eq. ()35j) . The total one-quark light-cone wave function is a sum of two contri- 
butions: 

F((p ± , z) = F[^(p ± , z) + F[^(p ± , z) (39) 



4 Distribution amplitudes of N and A 

Distribution amplitudes are, by definition, those components of the Fock vector of state 
for a given particle, which contain the lowest possible number of partons. For the nucleon 
the distribution amplitudes are its three-quark wave function on the light cone. As it is 
well-known, distribution amplitudes describe properties of the nucleon in hard exclusive 
processes. 

The full wave function of the nucleon is given by eq. (j!8|) and the wave function of 
the lowest component (with N c quarks) is the product of N c wave functions F[(p±,z) of 
eq. ()39|) . However it is not the end of the story. The point is that the minimum in eq. (fTo^) 
for the pion field, which corresponds to the nucleon, is degenerate. Any pion field which 
is the flavour (or space) rotation or the translation of the field of eq. (fTK|) also gives a 
minimum of the action. In other words, integration over the pion field in the functional 
integral has zero modes which should be taken into account exactly. 

It is well-known [T3*l IT%1 IT9*] that the integration over translational zero modes leads 
to conservation of the momenta. In the context of our calculation of nucleon distribution 
amplitudes, this leads to the condition z\ + . . . z Nc = 1 and pi± + . . -Pn c ± = 1- As to the 
integration over flavour rotations, it gives rise to the quantum numbers of the nucleon- 
soliton: the state with given quantum numbers (in SU(2) flavour group they are: spin 
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J, isospin T and their projections J3 and T 3 ) is obtained as a projection of the rotating 
soliton on the definite rotational wave function (which appears to be Wigner D-function). 6 
As a result, accounting for rotational and translational zero modes, we obtain the 
following expression for the nucleon wave function: 



• • • z Nc ,p 1± . ..p Nc± ) =Aj dRV2JTT(-l) T+T W^>(R) 
■d Zl d 2 Plx [dz Nc d* VNcL5 ^ ^ _ ^ (27r)2(5(2) ^ 



x 



(2tt) 2 "J (2vr 



,Pi 



x R f 9 \F 9 x l(p 1± , Zl )... R f g N N: F 9 x » : (p NcX , z Nc ) (40) 

Here *&j^t 3 * s the wave function for some state from the rotational band of the soliton, 
R is a rotational matrix from the flavour group (we consider flavour group SU(2) at the 
moment), A is the normalizing coefficient. 

The integration over rotational and translational modes breaks down the factorization 
of the nucleon wave function into the product of one-quark wave functions. This is to be 
expected as rotation and translation are, strictly speaking, specific correction in N c . Let 
us note that the wave function of eq. (J4U)) is symmetric under the exchange of quarks, as 
it should be, since the colour wave function of this state (omitted in this expression) is 
completely antisymmetric and, thus, corresponds to the colourless baryon. 

In the realistic applications one has to put N c = 3 and discuss states with J = T = 1/2 
(nucleon) and J = T = 3/2 (A- resonance). The full wave function of these states is given 
by eq. ([40)1 . However usually we are interested in the so-called distribution amplitudes 
which are the wave functions integrated over all p±: 

<P(z 1 ,z 2 ,z 3 ) = 5(z 1 + z 2 + z 3 -l)6(z 1 )9(z 2 )9(z 2 ) J . . . ^{z 1} . . .)■ (41) 

By definition, the normalization of the distribution amplitude is chosen in such a way 
that the total integral over all three Z{ is equal to unity. 

It is known [26, 2T] that relativistic invariance and the symmetry considerations re- 
strict the general form of the quark distribution for the nucleon in such a way that it 
depends on two scalar functions only: completely symmetric <f> a (zi, z%, £3) and antisym- 
metric 4> s {z\, z 2 , Z3). For example, for the proton with spin up we have: 

mz u z 2 ,zz) = ( ™?hp. z ^ (2\u t d I u t> -\u t u i d t> ^d t u | u |>) + 
v 6 

+ UZl ^ Z3) (\uUld1>-\dUlu t>) (42) 

Straitforward calculations based on eq. (j4U|) . of course, reproduce this general structure 
and give the concrete expressions for the symmetric and antisymmetric parts of the dis- 
tribution amplitude. However, these formulae are somewhat cumbersome and we will not 

6 See, for example, Review El f° r the discussion of the rotational quantization of baryon-soliton. 
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write them here. As an example, let us give only the expression for valence quark (level) 
conribution to the symmetric and antisymmetric parts of the quark distribution: 



4> s (zi, z 2 , z 3 ) = c J d x± [2g(z 1 ,x±)g(z 2 , x±)g(z 3 , x±) - 

-x]_j(z 2 , x ± ) {g(z 1 ,x ± )j(z 3 , x ± ) + g(z 3 , x±)j(zi,x±))] , (43) 

*.(*,*,«0 = ^/«P«xxlj(*,xx)(«(*,« J JJ(*«x)- 

-g{z3,x±)j(zi,x±)), (44) 

where g and j are the Fourier transforms of the level wave functions with I = and I = 1 
in the longitudinal direction: 



g(z, x±) = J dx 3 exp (-ix 3 (zM N - e)) 



h(r) - —j(r) 
r 



j(z,x±) = j dx 3 exp (— ix 3 (zM.N — e)) (45) 

and Co is a normalizing constant. 

In the non-relativistic limit, j(r) g(r) and only the symmetric wave function <p s (z) 
survives. The antisymmetric part is completely due to the relativistic effects and sea 
quarks. 

The total wave function should be calculated with the account for sea quarks according 
to eq. pS)l. However, substituting the level and pair wave functions into eq. (J3~%j) we arrive 
at the integral which diverges at small z' due to the phase volume factor Jz/z' . 

In fact, we have already met the singularity of this type in the calculations of pion 
wave function (221 and generalized distributions [2D]- The reason for this singularity is 
simple: at small z the quark virtuality 

2 M 2 +p l± 
Pi ~ (46) 

z% 

becomes large and we cannot neglect any longer the dependence of the dynamical quark 
mass M(p) on this virtuality. The dependence of the constituent mass on the virtuality is 
known in the instanton vacuum model JUE1E|; substituting M(p) in eq. (|3*H|) by this 
function we obtain the converging sea quarks contribution which tends to zero at z — » 0. 
As to the level contribution it remains non-zero at small z. 

The total wave function accounting for the sea quarks can also be represented in the 
form of eq. (J44j) where the functions g j receive contributions both from valence and sea 
quarks. Plots of the symmetric (f) s (zi, z 2 , z 3 ) and antisymmetric (f) a (zi, z 2 , z 3 ) parts of the 
nucleon wave functions are displayed in Fig.l and Fig.2, correspondingly. It is seen from 
these curves that the antisymmetric part of the wave function is two orders of magnitude 
smaller than the symmetric one. This fact is of a numerical origin, parametrically both 
parts of the wave function are of the same order. 

In fact, the quark distribution amplitudes calculated here are valid only in the region 
zN c ~ 1. It is the main region for the nucleon wave function. However, approximations 
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Figure 1: Symmetric part of the nucleoli quark distribution amplitude <f> s (zi, Z2, 23) as a 
function of z-i at different z\. 
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o 




Figure 2: Antisymmetric part of the nucleon quark distribution amplitude (f> a (zi, z 2 , £3) 
function of z 2 at different z\. 
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used in the present paper break at the end-points, at z% = and Zi = 1 but for different 
reasons. 

In the z — > 1 limit one quark is carrying almost all nucleon momentum and the 
momenta of all other quarks are necessary small. This limit contradicts the large N c - 
approximation which leads to the factorized wave function corresponding to independent 
quarks. The z — > 1 asymptotics corresponds to rather rare configuration of the pion field 
and its contribution to the quark distribution is exponentially small in N c . 7 

Wave functions calculated here are neither valid in the limit z — > 0. The main reason 
is that at small z the virtuality of nucleon constituents become large and we are driven out 
of the region of the applicability of chiral effective lagrangian of eq. (0). We have already 
mentioned that we have to take into account dependence of the dynamical mass of the 
quark M(p) on its virtuality. However it can be shown that this effect does not lead to 
nullification of the quark distributions at z = 0. It seems that more important is an effect 
of gluons which are known to be produced intensively at small z, at least in perturbation 
theory. As a consequence, the relative contribution of pure 3-quark component decreases 
at small z. As it was already mentioned, all effects of this kind were neglected here due 
to the instanton vacuum parameter (Mp) 2 <C 1. 

The small antisymmetric part of the nucleon wave function contradicts strongly to the 
commonly used parametrizations of quark distributions which are obtained on the basis 
of form factor data or QCD sum rules (see, e.g., |27J). Even more significantly our wave 
function differs from the Chernyak-Zhitnizky's one [2E1 which they suggested many 
years ago analyzing QCD sum rules. Their wave function has strong asymmetry z\ <-> Z2 
and even changes sign as a function of Z\ and z 2 (This looks strange for the wave function 
of the nucleon which is the ground state in the sector with nonzero baryon number). Our 
nucleon wave function is much closer to the so-called asymptotic wave function (valid at 
arbitray high normalization point) which is 

(j> s {z x , z 2 , z 3 ) = 120 z x z 2 z 3 , <Pa(zi> z 2, za) = 0; (47) 

but still differs from it. 

Data on the asymmetry of the nucleon wave function were reanalyzed recently in 
Ref. in the framework of the QCD sum rules on the light-cone. According to the 
results of this analysis, the nucleon wave function is not very far from the asymptotics. 
Let us also note that the data on 7r-meson photoproduction on the threshold [10J also 
favor symmetric wave function. 

Rotational symmetry of the nucleon-soliton allows one to calculate immediately quark 
distributions for A-resonance as well. For A + -resonance with J 3 = +1/2 we obtain 

$fl£=l/2(*l> 2& > z *) = 7^0a(2i, z 2 , z 3 ) (\u t u | d f> +\u | d I u T> - 

-\d]u[u]>) (48) 

7 In principle, this situation can be also described by semiclassical methods but one has to find a new 
minimum in the functional integral over pion field corresponding to the limit z —> 1. It can be shown 
that corrections to the mean pion field become essential when parameter l/iV c ln(l — z) > 1. Let us note 
that this situation is rather general and one faces the same problem also in the calculation of structure 
functions at large N c ^5] . 
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In other words A-resonance is characterized only by one symmetric wave function 
<Pa(z 1 ,z 2 ,z 3 ) = -j=J d 2 x± [g(z 1 ,x±)g(z 2 ,x±)g(z 3 ,x±) + 

+x 2 ± j(z 2 , x±) (g(zi,xx)j(z 3 , x±) + g(z 3 , x±)j(z 1 ,x±))^ (49) 

The main source of the experimental information on the nucleon wave functions is the 
asymptotics of form factors at large Q 2 . Four form factors are measured: magnetic form 
factors of proton and neutron, axial nucleon form factor and the transitional A — > N one. 
In the region of hard Q 2 , the process is factorized into a product of the hard QCD part 
and corresponding wave function (see, e.g. j^SJE] )• Asymptotics of the form factors are: 



9 P = Q Gm p = 2/ J[dx][dy] [T 1 (f) s (x)(j) s (y) + x <-> y] , 

g n = Q 4 G Mn = -y J[dx][dy] [(Ti - T 2 )<t) s (x)(j) s (y) + x «-> y\ , 

9a = Q 4 9a(Q 2 ) = y J[dx)[dy] [(AT, + T 2 )cf> s (x)<f> s (y) + x <- y] , 
2\/2f r 

g Ap = Q a G Mp a + = J[d 3 x] [d 3 y] [(T x - T 2 )<j> s {x)<j> A {y) + x ^ y] (50) 

(we give here the simplified version of hard kernels in the assumption that the antisym- 
metric part of the nucleon wave function is zero). Here [dx] denotes the integration over 
the fractions of total momentum with the condition x\ + x 2 + 23 = 1, / = (167ra s /9) 2 
and functions T% and T 2 are equal to 

r 1= 1 1 1 

2 3 (1 - 2i) 2 y 3 (l - yi) 2 x 2 {\ - xi) 2 y 2 {\ - yi) 2 

1 



2223(1 - x 3 )y 2 y 3 {l - y ± y 

T2 = 7Z \ 7, r- (51) 

2i2 3 (l - xi)yiy 3 {l - 1/3) 

Substituting our expressions for the wave functions into eq. (jHUj) we see that integrals 
over Zi are divergent due to the end-point singularities both in the region Zi — 1 and 
Zi — 0, where our approximations done are not valid (see above). As to the divergence at 
Zi — 1, it does not look serious: we know anyway (from asymptotics of structure functions) 
that wave functions should behave at least as (1 - 2) 3/2 . The singularities are weak and 
after any appropriate regularization the contribution of this region to the form factors is 
negligible. 

A discussion of the singularities at 2 = is much more involved. First, we do not know 
any strict QCD theorem which states that the wave function should be zero at z = 0. 
The divergence is linear at N c = 3 but it becomes worse at larger number of colours. As 
a result, the main contribution to the integral comes from this region. 

Thus, the considered form factors are rather bad objects to be treated in the framework 
of our model. We cannot calculate absolute value of the form factors but can say only 
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Ratio 


Ch.-Zh. [28J 


G.-St. [HU] 


Large N c 
(this paper) 


Experiment 


9n/9p 


-0.48 


-0.1 


-1/3 


-0.45 ±0.1 jSHl 


9 a/ 9 P 


1.53 


1.13 


4/3 


i.35 mu 


9 p a/ 9 P 


0.01 


0.81 


v/2/12 


< 0.3^21 



Table 1: Form factor ratios. 

that they should be amplified by a large parameter (connected in turn with parameter 
(Mp) _1 ). However we can calculate ratios of form factors using the fact that divergent 
part is universal for all form factors. In fact, this prediction is based only on the limit of 
N c —>■ oo and it does not use any concrete realization of the nucleon-soliton model. 

Ratios of the form factors in the large N c limit in comparison with results of the 
two popular parametrizations of the nucleon wave function are presented in Table 1. 
More extensive and complete analysis of the phenomenological implications of baryon 
distribution amplitudes can be found in |HSl EH ESj • In these works the corresponding 
analysis has been performed in the framework of so-called heterotic approach. 

We see that, qualitatively, large N c results agree with data 8 . Unfortunately the quality 
of the data is still rather bad. For recent comprehensive reviews of theory and phenomenol- 
ogy of hard exclusive reactions see Refs. [UHl EZj • 

5 Conclusions 

We have shown that the limit of large number of colours in QCD allows one to tell much 
about the light-cone wave function of the nucleon. In this limit the nucleon is a soliton 
and its quark wave functions are almost completely determined by the requirements of 
the large- N c factorization. One needs to know only the wave function of the discrete level 
in the self-consistent meson field and one-quark, two-quark, etc. wave functions of the 
mesons which form this mean field. Then the factorization, which is valid at N c — * oo, is 
enough to calculate the whole light-cone wave function of the nucleon and, in particular, 
its 3-quark component — the distribution amplitude. 

Put in such a way, this statement looks like a strict theorem. In the present paper we 
realized this program for the concrete model of the nucleon structure of eq. (0) motivated 
by QCD instanton vacuum. The main drawback of the model is, of course, the fact 
that the gluon components of the wave function are parametrically suppressed (by the 
parameter (Mp) 2 <C 1) while we are used to think that gluons play an essential role in 
the nucleon structure. 

One should take into account that the distribution amplitudes calculated in this paper 
are normalized at a very low normalization point (p < 600 MeV). Evolving these functions 
even to 1 GeV will lead to significant amount of gluons in the wave function. Of course, 
calculations of the gluon components of the wave functions (as corrections in the density 

8 Let us note that asymptotic wave function contradicts the data strongly and is not displayed here 
for this reason. 
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of instantons) are highly desirable. We plan to discuss the role of gluons in the nucleon- 
soliton wave function in the separate publication. For recent detailed works on evolution 
equations for baryon distribution amplitudes see Refs. jUJ H2] • 

However, we believe that even in the leading order, the large N c nucleon wave func- 
tion is interesting by itself. It was derived in the relativistically invariant field theory and 
therefore it obeys all general theorems and sum rules. It can be checked that the resulting 
wave function reproduces correctly the structure functions of the nucleon calculated in 
Ref. |18J and the generalized distributions of Ref. [213 • This is to be expected as they 
are calculated in the same model, under the same approximations. However, this demon- 
strates one more time that the model is self-consistent. On the other hand it is known 
that the structure functions obtained in Ref. [TH] describe the data, at least qualitatively. 

The main features of the nucleon wave functions considered in this paper are: i) 
relatively small contribution of the sea quarks to the component with the lowest number 
of partons ii) soft distribution in the transverse momenta of quarks with average value 
« 400 MeV iii) the form which is much closer to the asymptotic wave function than to 
the wave functions proposed on the basis of QCD sum rules j2H] ( n o zeroes or change 
of a sign) but which is still rather far from the asymptotic form, iv) wave function is 
singular (does not vanish) at the end-points. These singularities can disappear as a result 
of evolution or if one takes into account corrections in N c . Nevertheless it is natural to 
expect that factorization in the baryon channel will appear at higher Q 2 than in meson 
sector. It is quite probable that this is really the case in nature [39J. 

Let us note that the wave functions of mesons in the instanton vacuum appears to be 
rather close to asymptotic one in agreement with the data (see Ref. [22] )• Baryon wave 
function in the same model is far from the asymptotics and this fact again corresponds 
to the data (asymptotics of the form factors). This is natural, as in the limit N c — > oo 
of the nature of mesons and baryons is completely different. We conclude that this limit 
describes experimental wave functions at least on the qualitative level. 
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